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We develop a flexible quasiparticle theory of transport coefficients of hot hadronic 
matter at finite baryon density. We begin with a hadronic quasiparticle model which 
includes a scalar and a vector mean field. Quasiparticle energies and the mean 
fields depend on temperature and baryon chemical potential. Starting with the 
quasiparticle dispersion relation, we derive the Boltzmann equation and use the 
Chapman-Enskog expansion to derive formulas for the shear and bulk viscosities 
and thermal conductivity. We obtain both relaxation time approximation formulas 
and more general integral equations. Throughout the work, we explicitly enforce 
the Landau-Lifshitz conditions of fit and ensure the theory is thermodynamically 
self-consistent. The derived formulas should be useful for predicting the transport 
coefficients of the hadronic phase of matter produced in heavy-ion collisions at the 
Relativistic Heavy Ion Collider (RHIC) and at other accelerators. 


INTRODUCTION 


A central challenge in nuclear physics is elucidating the structure of the Quantum Chro¬ 
modynamics (QCD) phase diagram. Based on theoretical models, it is widely-believed that 
the phase diagram contains a line of hrst-order phase transition which ends at a point of 
second-order phase transition - the critical point Despite a dedicated search for the 

critical point with the first beam energy scan at the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Lab, and possible hints of the critical point|3|, the location of the 
critical point remains a mystery. From lattice QCD calculations we know the transi¬ 

tion from hadrons to quarks and gluons is an analytic crossover near temperature T 150 
MeV at zero baryon chemical potential Hence, the critical point is generally thought to 
be located at T < 160 MeV and hb equal to several hundreds of MeV. 

A second, future beam energy scan at RHIC will search for the critical point with 
greatly increased statistics and upgraded detectors B. To maximize the discovery potential, 
experimental efforts must be accompanied by complementary improvements in theoretical 
modeling of QCD matter at moderate temperatures and large baryon chemical potentials. 
In p revious papers, we investigated the equation of state at finite baryon chemical potential 


m 11| . In this work, we derive new formulas to compute the shear and bulk viscosities 


and thermal conductivity of hot hadronic matter with /i^ > 0. We employ a flexible, 
thermodynamically consistent framework of hadronic quasiparticles with medium-dependent 
quasiparticle masses and with a scalar and vector mean field. This may be considered a 
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12| to include nonzero baryon chemical potential and the concomitant 


natural extension of 
vector mean held. 

Transport coefficients like the shear and bulk viscosities and thermal conductivity are 
especially interesting quantities to study for several reasons. First, they are essential theoret¬ 
ical inputs for hydrodynamic simulations, which are critical tools for interpreting heavy-ion 
collision data. In hydrodynamic simulations, the shear and bulk viscosities inhuence various 
observables, such as the elliptic how coefficients Vn and the hadron transverse momentum 
{Pt) spectrum [l^ 16|. Furthermore, the temperature and chemical potential dependence of 
transport coefficients may reveal the location of phase transitions: in many physical systems, 
the shear viscosity is a minimum and the bulk viscosity a maximum at the phase transition 
17l |. A third motivation is investigating the KSS lower bound 18| on the shear viscosity to 


entropy density t]/ s > l/dvr for strongly-coupled conformal theories and its implications for 
QCD. 

In principle, the transport coefficients can be computed directly from QCD using the 
Kubo formulas [l^. However, QCD is strongly coupled at energies accessible to heavy- 
ion collision experiments, complicating hrst-principles calculations. There were some early 
attempts to employ lattice QCD jl^, 21], but even today it is challenging to achieve a large 
enough grid with a small enough grid spacing to accurately compute transport coefficients. 
Furthermore, lattice QCD simulations are currently very difficult at hnite baryon chemical 
potential due to the well-known fermion sign problem. Hence, many of the early works 


22h 25| computed transport coefficients of quark-gluon plasmas, or hadronic gases with a 


few species of particles, using the Boltzmann equation in the relaxation time approximation. 
These early works did not include mean helds or medium-dependent masses. 

Later on. Jeon [ 2 ^ and Jeon and Yaffe computed the shear and bulk viscosities of a 
hot, weakly-coupled scalar held theory using perturbation theory. Amazingly, they showed 
that their complicated perturbative calculation of transport coefficients was reproduced by 
a simpler kinetic theory of quasiparticles with temperature-dependent masses and a scalar 
mean held. The same conclusion was found for hot, weakly-coupled QCD and QED 


This was also consistent with an earlier analysis of transport in a nucleon plus a meson 
system, which similarly found that renormalized quasiparticle masses were required js^. 
Though astounding, this makes intuitive sense: kinetic theory is widely used to model non¬ 
equilibrium systems, and renormalized particle masses are ubiquitous in hnite-temperature 
held theories. (They are also present in Fermi liquid theory [37|.) Also, temperature- and 
chemical potential-dependent masses allow quasiparticle models to generate more realistic. 


non-ideal gas, equations of state |38|. Furthermore, as Gorenstein and Yang pointed out 
3^, the scalar mean held is essential for maintaining thermodynamic self consistency when 
masses depend on temperature and/or chemical potential. Hence, it seems kinetic theories 
of quasiparticles with medium-dependent masses and mean helds are powerful theoretical 
tools, though thermodynamic consistency must be carefully maintained. 

More recentl y, t he conjecture of a lower bound on rj/s by Kovtun, Son, and Starinets 


from AdS/CFT M 
calculations l40l-l43l 


ignited a hurry of additional work. There were several more lattice 
There were also many studies with Boltzmann equations - most of 
them without medium-dependent masses or mean helds. Shear viscosity was computed for 
pion-nucleon gases at low temperatures and varying chemical potentials in 44. l45l . Bulk 
viscosity of cool pion gases was computed using chiral perturbation theory in 45 J46| . Shear 
viscosity in mixtures of hadrons with excluded volumes were calculated in [47-49 . 

There were a few attempts to employ the more powerful quasiparticle models with 
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medium-dependent masses to compute transport coefficients. In an early work, Sasaki and 
Redlich applied kinetic theory and the relaxation time approximation to a quasiparticle 
model to compute the bulk viscosity near a chiral phase transition [s^. Later, Chakraborty 
and one of us developed a comprehensive theory of shear and bulk viscosities in hadronic 


gases [12|. That work included multiple hadron species with temperature-dependent masses 


and a scalar mean field in a thermodynamically self consistent way. They derived formulas 
for shear and bulk viscosity and provided both relaxation time approximation formulas and 
more general integral equations. However, they did not include chemical potentials, hence, 
thermal conductivity was not considered in that work. Bluhm, Kampfer, and Redlich used 
a similar quasiparticle formalism to study the shear and bulk viscosity of gluon matter in 
511 1 (also without chemical potentials). Thus, a natural question is, how does the formalism 


of HI generalize to finite baryon chemical potential? Also, what is the formula for thermal 
conductivity? 

Several papers have tried different ansatzes for generalized viscosity formulas (in the 
relaxation time approximation) when the baryon chemical potential is non-zero. Chen, 
Liu, Song, and Wang calculated the shear and bulk viscosities of weakly-coupled quark 
gluon plasma at finite temperature and chemical potential in j52| using a quasiparticle 
model with medium-dependent masses and a scalar mean field. Khvorostukhin, Toneev, 
and Voskresensky compared three ansatzes for the generalized bulk viscosity formula 15311 
of a hadron gas with medium-dependent masses and a scalar mean field; see also 5J, l55l . 
Interestingly, Khvorostukhin’s quasiparticle model also included a vector (a;) mean field [53l. 


54| : as is well known, they are important to account for repulsive forces in hadronic matter 


with large baryon densities. This type of model is quite relevant for studying the moderate 
temperature hadronic matter formed in the beam energy scan at RHIC. It is also relevant for 
experiments at the Super Proton Synchrotron Heavy Ion and Neutrino Experiment (SHINE) 
at CERN and at the future Facility for Antiproton and Ion Research (FAIR) at GSI. Given 
the usefulness of this kind of model, it is desirable to put the results on a firmer theoretical 
foundation and (ideally) determine which of the ansatzes presented in [HI and 531 are 
correct. 

In this work, we present detailed derivations of the formulas for the shear and bulk 
viscosities and thermal conductivity of a gas of hadronic quasiparticles. We include a scalar 
and a vector mean field, where the mean fields and the quasiparticle masses depend on 
temperature and baryon chemical potential. Generalization to multiple scalar and vector 
fields is straightforward but not included here for clarity of presentation. Starting from 
the quasiparticle dispersion relation, we obtain the Boltzmann equation, and then use the 
Ghapman-Enskog expansion to derive formulas for the transport coefficients. At each step 
we ensure that thermodynamic self-consistency is maintained, and we carefully enforce the 
Landau-Lifshitz conditions of fit; we later show this is vital to obtaining the correct results. 
We derive both relaxation time approximation formulas and more general integral equa¬ 
tions. Finally, we show that the formulas for shear and bulk viscosities are straightforward 
generalizations of previous results [HI, [13 if one recalls that entropy per baryon is conserved 
in ideal hydrodynamics (neglecting viscous effects). Glassical statistics are used in the main 
text for ease of presentation, but results which include quantum statistics are presented in 
the appendix, albeit without detailed derivations. 
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II. QUASIPARTICLES 


In this section we discuss quasiparticle dispersion relations for baryons and mesons. In 
the simplest mean held approach all hadrons acquire effective masses in the medium. In 
addition, baryons acquire effective chemical potentials. We will focus attention on baryons 
since the inclusion of the baryon chemical potential is the new feature of this work compared 


to 12 


The piece of the Lagrangian involving baryons is 


£ 


baryon 




rrij + - g^j 


( 1 ) 


Here j refers to the species of baryon. For simplicity of presentation we include only a generic 
scalar meson a and a generic vector meson u. When evaluating the partition function there 
enters an additional term of the form where /tb is the baryon chemical potential. 

Since we are using Dirac spinors both particles and antiparticles are included. Particles have 
chemical potential hb while antiparticles have chemical potential —^b- 

For a uniform medium in thermal equilibrium the meson fields acquire space-time inde¬ 
pendent nonzero mean values denoted by a and in the rest frame of the medium the 
spatial part of the vector field vanishes on account of rotational symmetry, ch = 0, but in a 
general frame of reference it does not. The dispersion relation for particles is 


E+ (p) = y/(p - gu^jU^y + mf + g^juj^ (2) 

and for antiparticles 

EJ (p) = y^(p + gujj^Y + mf - g^jU^ . (3) 

The kinetic momentum p* is related to the canonical momentum p by 

P*j = P- (4) 

for particles and by 

Pi = P + (5) 

for antiparticles. Particles and antiparticles have a common mass m*. In this mean field 
approach it is given by m* = rrij — g^jd'. 

A more convenient way to think about the dispersion relations is to recognize a shift in 
both the mass and chemical potential of quasiparticles and anti-quasiparticles. They both 
have energy 

E*^{p*) = p*^ + mf (6) 

while their chemical potentials are opposite in sign 

(7) 

as behts particles and antiparticles. 

Mesons do not have a baryon chemical potential. They could have chemical potentials 
for electric charge or strangeness, but we do not consider that possibility here for simplicity. 
Hence their dispersion relations, in mean held approximation, are of the form 

E*{p) = \/p^ + . 


( 8 ) 
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Note that the kinetic and canonical momenta are the same for mesons. The effective masses 
and effective chemical potentials can be found self-consistently once one hxes the Lagrangian. 
In equilibrium the phase space density for a particle (or antiparticle) of type a is given 


by 


/a(x, p*,t) 


1 


(9) 


Here Sq denotes the spin. There are Fermi-Dirac and Bose-Einstein distributions. Later 
on we will simplify our results by using classical statistics, although that approximation is 
not necessary. Results including quantum statistics are given in the appendix. Momentum 
space integration will be abbreviated as 


dr 


* 

a 


{2Sa + 1) 


dYa 

(27r)3 


( 10 ) 


indicating that the kinetic momentum is chosen as the independent variable, and the spin 
degeneracy is included. 


III. BOLTZMANN EQUATION 

The general form of the Boltzmann equation for the distribution function /^(x, p*,t) is 


dt ’ ^ ’ dt dx^ dt dp*^ dt 


= Cn. 


( 11 ) 


The right hand side is the collision term which will be discussed later. Here we focus on the 
left hand side. It involves the trajectory x(f) and p*{t) between collisions. This trajectory 
is in general not a straight line because the particle is moving in a mean held which can be 
space and time dependent. 

The velocity is 

dx* dEa p” ^^2) 


dt dp\ El 

The relativistic version of Newton’s Second Law is 


dt 




Note that it is p that is held hxed, not p*. The right hand side is 


dx^ 


ml dml 
E* dx^ 


9i^ 


‘ dx^ E* 


+ 9o. 


du° 
‘ dx^ 


(13) 


(14) 


The left hand side of Newton’s Second Law can be written in terms of the kinetic momentum 
as 

Pn* dp*^ du}'- dp*^ /9a}* p*^ dui^\ 

(15) 


dpi 

dt 


dt 


** dch* _ dp** 

T 9uja TV 77 1“ 9uia 


dt 


dt 


duj'' p*l 9a;* 

+ 


dt El dx^ 


The time derivatives of x and p* can now be replaced in Eq. (ITTl) to put the Boltzmann 
equation in the form 

dfa , * dfa p** dfa 

-(^,p ,t) = — + - ^ 

















6 


dfa f m* dm* , ^ 


dcj^ du^ p*^ f du^ du^ 

dx^ dt ^ i?* V dx^ dx^ 


dp*^ I E* dx^ 

This can be simplified by making use of the kinetic 4-momentum 


= Cn 


and the field strength tensor 


Pa 




The final form is 


dfa, * . P*^ ^ r 

^(x.p .t) = 


E* dx^ ^^'"E* 


dfa 

dp*i 


Ca. 


(16) 

(17) 

(18) 

(19) 


IV. ENERGY-MOMENTUM TENSOR AND BARYON CURRENT 


In this section we present the structure of the energy-momentum tensor and of the 
baryon current In terms of temperature, chemical potential, and flow velocity they 
are 


= -Pgi^'^ + wu^^u'^ + 


( 20 ) 


and 


= hbu^ + AJ ^ ( 21 ) 

where P{T,pb) is the pressure, s = dP/dT is the entropy density, ub = dP/dpB is the 
baryon density, e = —P + Ts + pb^b is the energy density, and w = e + P \s the en¬ 
thalpy density. In the Landau-Lifshitz approach, which we use, is the velocity of energy 
transport. The AT^^ and AJ)^ are dissipative parts given by 


and 


AT^" = p + fA'^^aX) - 


AJg = A X) (f) . 


M'' 


( 22 ) 


(23) 


Here p, ( and A are the shear viscosity, bulk viscosity, and thermal conductivity, respectively. 
The other symbols are 

D = uPdp , (24) 

= d^^ - u>^D , (25) 

. (26) 

Our metric is (-f, —, —, —). Additionally, the entropy current is 


3“ = su^ - J" . 


(27) 


Now we need to express T'^'' and J)) in terms of the quasiparticles and mean fields. One 
expression for the former is 


T'"' = J2 j 


* Fa Fa 


E* 


/„ + t/""!/ (<T, Fu]f) + mlu}V. 


(28) 
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The first term is familiar as the kinetic contribution. The second term is the usual meson 
field potential energy; it includes the mass terms and — plus any interaction 

terms which are more than two powers of the fields. Note that kinetic terms for the mean 
meson helds are not included because they are second order in space-time gradients and are 
not included in first order viscous fluid dynamics. The last term is not obviously of the form 
of Eq. 0201) . However, when one remembers that is the energy flux in the direction i, 
and that Ea is the complete quasiparticle energy and not E*, then one would write 

r"" = E / + 3‘“'U{iT, w'-w,). (29) 


Using pi^ = p*J^ we get 

+ g^^U{a, oo^oop) 9-- / ^r:|y/a • (30) 

a ^ a ^ 

The vector mean field is determined by its equation of motion. Assuming an interaction 
only with the baryons (this assumption is easily relaxed) it is 

{d^ + ml) ^ (31) 

3 

where the averaging refers to the quasiparticle distribution. Recognizing that the summation 
index j refers to both baryons and antibaryons, and dropping the d’Alembertian because of 
hrst order viscous fluid dynamics, we have 

= E 9“(32) 

n ^ 


(We remind the reader that the coupling g^^a is opposite in sign for baryons and antibaryons.) 
Hence Eqs. (128|) and (|^ are the same. 

In a similar way the scalar mean field is determined by its equation of motion. This 
turns out to be 


ad 

a 




dK^fc 

^ E* ^ 


(33) 


The coupling to scalar mesons of baryons and antibaryons has the same sign, unlike the 
coupling to vector mesons. 

The structure of the baryon current is readily deduced to be 



(34) 


where ba denotes the baryon number of a. 

It can be shown that energy and momentum are conserved, namely 


dpT^’' = 0 


(35) 


9„Jg = 0. 


and so is baryon number 


(36) 





These conservation laws follow from the requirement that 


f rfr;x.c<. = 0, (37) 

a 

The Xa represents the contribution from quasiparticle a to any conserved quantity, such 
as energy, momentum, or baryon number. The calculations are straightforward but very 
lengthy and tedious. We have performed them but they are not reproduced here. It is also 
straightforward, and much less tedious, to show that the mean held equation of state follows 
from the above expressions for and Jg when the system is uniform, time independent, 
and in thermal and chemical equilibrium. 


V. DEPARTURES FROM EQUILIBRIUM OF THE QUASIPARTICLE 

DISTRIBUTION FUNCTION 

To hrst order in departures from equilibrium, we can express the quasiparticle distribu¬ 
tion function as 

h = (1 + 0a) (38) 

where is the distribution function in thermal and chemical equilibrium. The nonequilib- 
rium part (j)a leads to the nonequilibrium contributions and A J^, so 0a must contain 

the same space-time gradients as found in them. Therefore, 0a must have the form 

0a = -AadpU'' - BaPaDy + ^aPaPa + D^Up + ^Ap^dpUP) . (39) 

The functions Aa, Ba and Ca only depend on momentum p while only depends on space- 
time coordinate x. 

The departure from equilibrium of the quasiparticle distributions can be used to compute 
the departure from equilibrium of the energy-momentum tensor. It is convenient to work in 
the local rest frame. The variation of the space-space part of expression fl28l) is 

I ■5/. - ) + g'’su. (40) 

u. d \ d / 

To obtain the variation in the mean field potential we start with the expression for the 
pressure P{T,pb) = Bq — U. Here Pq is the kinetic contribution to the pressure from the 
quasiparticles. The entropy density is obtained from s = dP(T, pB)/dT. This has three 
contributions: the first is from Sq which is the same functional form as for particles with T- 
and /ifi-independent energies, the second is from the variation of the quasiparticle energies 
due to variations in T and ps, and hnally there is the contribution —dU/dT at fixed pb- 
The mean field carries no entropy, therefore the second and third terms must cancel. Using 
classical statistics for simplicity we have 

p„ = t 5^/dr:/:- (41) 


E 


/ 


dr: 


dEA 

dT ) 


Ms 


and thus 


dU 


req 

J d 


(42) 
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The same argument applies to differentiation with respect to ps, which gives the baryon 
density. The mean held carries no baryon number, so similarly 


dU 

d^B 


E 


dr: 


d^B ) T 


(43) 


Hence 

iU = -Y, f dr'JEaf^ (44) 

a 

where Ea = and 

TH * 

SEa = ■ (45) 

E*a 

Now we come to the deviation in the quasiparticle distribution function. The fa in 
general will have departures from the equilibrium form, but it can also change because the 
quasiparticle energy departs from its eqnilibrium value. Let us denote E^ the eqnilibrium 
valne and Ea the total noneqnilibrium energy; it is the latter which is conserved in the 
particle collisions. Similarly, we denote T° and /i^ the eqnilibrium values. Then we write 


= fT(Efi,T^,A) + 6U, 

/.(£■,, T, ij-b) = C(Ea, T°, 4) + Sfa. (46) 

The deviations are related to each other by 

4/a = 4/a + (1^) ^ JE, = Si - (47) 

where the second eqnality follows when using classical statistics. 

It is always the 6fa which determine the transport coefficients. Therefore we express 
(5T®-^ in terms of 6fa instead of 5fa- 


= 


Y. I “K^Wsl -YEr 


El 


4 Ira yg 




E* 


/•eq 
J a 


E’Y / 




(48) 


The integrand of the second term depends only on the magnitude of p*, apart from the 
factor Therefore, one may effectively make the replacement p*fp*J —)■ Then 

the terms not involving 6fa all have a factor of They can be written as a sum of 


~T~ 



dr* T- 


IP. 


* 12 


3e: 


re 
J a 


eq 


and 




6m*„ 




IP 


* 12 


Ip 


* 12 


3TE* 


3Ef 


/, 


eq 
a ‘ 
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It can be shown that both of these integrate to zero (using classical statistics). Hence we 
find 

(49) 

a “ 

as our final result. 

The variation in the time-time component of the energy-momentum tensor, starting with 
either Eq. (1281) or (|29|) . is 

[ dTlEJfa. (50) 

a 

We use Eq. fITTD for 5fa- The variation of the local energy Ea 

^777*2 

5Ea = ^ + g^adu^ (51) 

2E* 


can be expressed in terms of the variations in temperature and chemical potential 


5mf 

5u^ 




dT ) 


5T 


Ms 


df^B A 


5jlB ■ 


(52) 

(53) 


The variations 6T and djjLB are not independent. They are related by the hydrodynamic 
flow of the matter which to this order occurs at constant entropy per baryon a = s/ub- 
Dissipation should not be included since it would lead to second-order effects which are 
consistently neglected in first order viscous fluid dynamics. The relation can be expressed 
in various ways, including these: 


Here vf 


1 / y 

\ifr ), Tv'i tL 

{dP/de)x is the speed of sound at constant x. It 


_ Xtt — crXiiT 

(yX^B^l - XfiT ' 

is easily shown that 


(54) 


2 _ ^BXfjT 

T{xttX/b)i - xIt) ’ 

^2 _ ^BXTT — SX^iT 

hs(XTrX/2At “ x^r) ’ 
2 vITs + vI^iub 


(55) 


relationships that are independent of the specific equation of state. Of course waves do 
not physically propagate at constant n or s, only at constant a, but these definitions are 
useful for various intermediate steps in various applications. The other symbol represents 
the susceptibilities 


Xxy 


d^P(T,B) 

dxdy 


(56) 
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Rather than thinking of m* and a;® as functions of T and /is we can think of them as 
functions of T and a. Then 

(¥) "■ 

,T. 


Smf = 


= 


dTj^ 


(57) 

(58) 


Next, we need to relate the variations in T and /is to the variation 6fa- The latter 
variation is done at hxed En and is 


5 fa = /, 


eq 


Ea — + T 


dfia\ 

dT I 


J'2 


(59) 


(Recall that Ha = baHs-) The term from Eq. (1471) which needs to be rewritten is 

■T2 {dmf/dT^)^ + g^aT {du^dT)^ E^ 


rji J a 


Ea- Ha + T {dfia/dT)^ 
T{dEa/dT), 


Sfa 


Sfa 


(60) 


[Ea-Ha + T{dHa/dT)^ 

We reiterate that the temperature derivative of a function E depending on T and /is, taken 
at fixed entropy per baryon, is 


dE\ _ fdE\ 


dT \dT 

/ <J \ / i 


MB 


dfiB/j 


dT J \dTj 

y <7 N / , 


Ms 


T v‘^ 


<9/is Jj 


The hnal expression is therefore 

dTlE, 


T{dEa/dT), 


Ea-Ha + T{dHa/dT)^ 


5 fa. 


(61) 


(62) 


When the baryon density goes to zero this reduces to the formula known in the literature. 
The time-space component has the very natural form 


„ J a 


(63) 


To express this in terms of 6fa, we note that the last term on the right-hand side of Eq. 
(ITT)) is spherically symmetric in momentum space and therefore that term integrates to zero. 
This is not true of the other term because the deviation (pa does have terms that depend on 
the direction of the momentum. Therefore 


ST«’ = yI . 

n Q. 


(64) 


Lastly we need the variations in the baryon current. The steps are by now very familiar. 
The results are 

T{dEa/dT)^ 


5Jl = 




dv: 


Ea- Ha+ T {dua/dT)^ 


5fa 


and 


sj'b = E 




( 66 ) 


(66) 
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VI. GENERAL FORMULAS FOR THE TRANSPORT COEFFICIENTS 


Suppose that we know the scalars Aa, Ba, and Ca in Eq. fl39l) as functions of the magni¬ 
tude of the momentum p*. Then in the local rest frame we should equate the hydrodynamic 
expression AT*-^ from Eq. (j22ll with the quasiparticle expression (fT*-’ from Eq. (091), the 
latter being 




E/- 


r 

* Fa Fa 


E* 


- AadaUP - BaP^D,. I 


T ) 


CaPaPa {DfaU^ + D^Uf, + ^Af^^dpUP) 


/, 


eq 


(67) 


The Ba integrates to zero by symmetry. In the local rest frame the derivative dkUo = 0, so 
the the summation over p and z/ is a sum over spatial indices kl only. In the Aa term we 
can use 

Pa Pa ^ 3 IPal "U 

and in the Ca term we can use 


P*>*aP*aPa 


-)■ 


\V‘*a\^{^ij^kl + Ak^jl + ^iAjk) 


because in the local rest frame p = p*. Equating the tensorial structures then gives ns the 
shear viscosity 

= ( 68 ) 

a '' “ 

and the bulk viscosity 

«= 5 E / . (69) 

a '' “ 

For the baryon current we compare the AJ^ from Eq. fl2^ with the dissipative part of 
Eq. fIMl) in the local rest frame. The latter is 


a 


dr 


'E* 


-BaP^D 



f 

J a 


eq 


(70) 


Obviously the Aa and Ca terms integrate to zero on account of symmetry. After some 
manipulation this results in an expression for the thermal conductivity 


A = 


1 

3 



(71) 


To solve for the functions Ba, and Ca we turn to the Chapman-Enskog method. 
This entails expanding both sides of the Boltzmann equation flT^ to first order in the (j)a- 
It leads to integral equations which in general must be solved numerically. 
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Here we follow the notation of [12(]. Including 2-to-2, 2-to-l and l-to-2 processes, and 
using classical statistics (these restrictions are easily relaxed) the collision integral is 


bed 


dri dr: <ir’jW(a, 6|c, d){fM - fM 


f dr: dr,W {a\c, d){fj, - f^} 

_ 7 


cd 


be 


The W are given as 

(2 

W (a, b\c, d) = — 

and 

W{a\c,d) = 


2E*2E*2E*2E* 


^4^4/ 


2E*JE*2E* 


fafb} ■ 

(72) 

Ad(a, b\c, d)p 

(73) 

-\M{a\c, d) ^. 

(74) 


The use of E* instead of Ea in the denominators ensures that the phase space integration 
is Lorentz covariant. Also note that, following Larionov 5^, we use dimensionless matrix 
elements Ai averaged over spin in both initial and hnal states. This is necessary to balance 
the degeneracy factors in the dT*. We use chemical equilibrium (for example, a + h c + d 
gives Then the collision integral becomes 


bed 


dV; dTl dT* CW{a, b\c, d) [0, + ct)d-^a- H 


+ ra^Yl / c?r:dr*w(a|c,d)[0, + 0,-0,] 

cd 

+ /r E / ‘'r; -ir: fpw {c\a, t) 14,, -4,- 4,]. 


(76) 


be 


This constitutes the right-hand-side of the Boltzmann equation. 

The left-hand side of the Boltzmann equation flTU]) is computed using the local equilib¬ 
rium form of the distribution function 


/:‘i(a:,p*)=exp 


T{x) 


exp 


T(x) 


exp 


'Xc(x)pT 

T(x) 


exp 


P'a(x) 

T{x) 


(76) 


Here the flow velocity, temperature and chemical potential all depend on x. Although not 
explicitly indicated, depends on x via the dependence of m* and a)" on x, while E* 
depends on x via m* only. The left-hand side must be expressed in terms of the same 
space-time gradients as (j)a, namely dpU^, Dy (fis/T), and [D^Uy + DyU^ + The 

calculation is long and tedious. Space-time derivatives of T and fiB are expressed in terms of 
the relevant tensor structures by using the perfect fluid equations for conservation of energy, 
momentum and baryon number. Some useful intermediate results are 


DT = -vlTdpUp, 

D^b = -v'^fiBdpUp. 


(77) 
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One form of the left-hand side (in the local rest frame) is 


> a _ £eq 

J a 

+ ( k 

-1 


dt 


Ip:p 

3TE* 


dy 


eg PaPa 


2TE* 


, ( Ea ha 

^ dT 


E ^Vju E Ei/U^ ^ 






(78) 


Now Ea — ha in the hrst line could be replaced by E* — p*, and Ea in the second line could 
be replaced by E* -|- g^ad}^- With a little manipulation this can be shown to be equivalent 
to Sasaki and Redlich who, however, did not include a vector held nor the D^^hb/T) term. 
Another form is to write out the derivatives in the hrst line explicitly. This results in 


dt 


_ /eq 
J t 


3TE* 


PaP 


E 


*2 _j2( dmf \ ^ 


V aT2 y 

N / ( 


dT \ T 


E* 


dpUp 




pK 


/r ^ ( Dpu, + D,u, + -^p^dpu 


(79) 


In the limit that the chemical potential goes to zero this reproduces the results of Jeon and 
Yahe 27| and of Chakraborty and Kapusta fl^ . 

Now we subtract the right-hand side from the left-hand side and set the resulting ex¬ 
pression to zero. This leads to 


A {dpup) + 


Pb \ 
T ) 


Cl” (Dpu„ + = 0 


(80) 


where 


A = 


3TE[ 

Y 

bed 


■k* 12 


3vt 


E. 


*2 


T2 


dm 


*2 


dT^ 


.2 d ffi: 




1 + (5, 


ab 


dTl dT*a dr: /r W(a, 6|c, d) A + Ad-Aa- A,] 


+ J2[ ^^r:dr:w(a|c,d) A-fArf-A,] 

cd d 

+ ^ / dr; dr; 4”’ w'(c|o. h) [A, -a^- a^] 


(81) 


be 























15 


and 


Bit = b. 


nBEa\ V^a 


w E, 


E 

bed 


1 6, 


ab 


dr; dr; dr; fp ^Ic. d) |b=p? + B^p; - b^p; - srfi 


/ dr: <ir: w(aic, d) |b,p: + BjpS - b.p:i 

cd 

Y, f dr; dr; fp w{c\a, b) |s,p: - b„p: - Sjp^] 


be 


and 


^ PaPa 
a 2E*T 


E 

bed 


(82) 


1 + (5, 


ab 


dr; dr; dr; c vr(o, 6|c, d) [c.pK + c^p^p; - c^ApI - CjpK] 


5^ / dr: dr: vr(a|c, d) [c.pJp: + c^pSp: - CM] 

cd 

Y f dr; dr; c w(c\a, b) [cm + cm - cm] ■ 


(83) 


be 


Dne to the tensorial strnctnre of these eqnations the solntion requires that = 0, = 0, 

and = 0. These are integral equations for the functions Aa, Ba, and Ca which depend 
on the magnitude of the momentum p*. 


VII. LANDAU-LIFSHTIZ CONDITIONS OF FIT 


The set of equations fl8T|) - fl83|) are integral equations for the functions Aa, Ba, and 
Ca- Consider the equation for Aa- If we have a particular solution A'^^^ we can generate 
another solution Aa = A^^ — aEEa — asba where the constant coefficients aE and are 
independent of particle type a. The reason is that energy and baryon number are conserved 
in the collision, decay, and fusion processes. This arbitrariness exists because of the freedom 
to define the local rest frame or, equivalently, the flow velocity u^. This arbitrariness is 
removed once one specifies the Landau-Lifshitz definition of the local rest frame, also called 
the condition of fit. Requiring that = 0 in the local rest frame results in 


ttE 




T{dEaldT)a 

^ia + T{dlialdT)^_ 


+aB ^ba f dT*aE, 
a 

= E / 


1 - 


T{dEa/dT), 


Ea- fia + T{dfia/dT) 

T{dEa/dT)a 


f 

J ( 


eq 


Ea-fIa + T(dfda/dT)^ 


4 par Teq 
■^a / a ' 


(84) 
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Requiring that (5Jg = 0 in the local rest frame results in 

T{dEa/dT), 


f 

a 

f 

a 


dVlE, 


1 - 


dr: 


dv: 


Let us express these equations as 


Ea-^Ia + T{^^ia/^T)^ 

TjdEjdT), 

Ea- IJ.a + T (dHa/dT)^ 

T{dEjdT), 


req 
J a 

feq 

J a 


Ea-fia + T{dfia/dT)J 


4 par fe 

■^a J a 


eq 


(85) 


O'eXe + — Ze , 

oeRe + ueLe = Zb ■ 


( 86 ) 


The solutions are 


Ye^e - XeZ 


Qb — 


Be — 


E^B 


YeXb - XeYb ’ 
XbZb — YbZe 


YeXb - XeYb 

When these are substituted into the expression fl6^ for the bulk viscosity we get 

1 


C = 3 E / dT*}^CA^ 

„ J a 


TubCLb — Twcle ■ 


(87) 


( 88 ) 


First consider the case where there are no mean fields, only on-shell particles traveling 
in vacuum and undergoing localized collisions. In this case 6fa = Sfa, and one finds 

Xe = T{T‘^Xtt + “^fiBEXfiT + ^^BX^l^l), 

Xb = T{TXiiT + 5 

Xe = T{Txht + hfiXMAt) • 

Xb = Txtat, (89) 

The combination of oe and qb which is needed for the bulk viscosity is 


TnBttB + Twbe = v^Ze + (n^ - vDubZ, 


B 


E / ‘'r; [vlE, + X: - 


(90) 


Here Ea = E* = ^/p'^ + because of the absence of mean helds. The bulk viscosity is then 
1 


C = 5E/< 

a 

This is a limiting form of 


Ip:p 

e: 


- 3 [vIE: + - vl)baPB] An 


ceq 


(91) 


Ip:p 

E* 


+3vn ^ 


dT 


T 


4 par req 
■^a J a ? 


(92) 
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once one recognizes Eq. fl5^ . This makes perfect sense because the modihcation of the 
integrand compared to Eq. (I6^ matches the structure of the source of Aa in Eq. (IHTj) . 

It is not easy to find simple expressions for X^, Xb, Ye, yn when mean helds are included, 
hence there are no simple expressions for and a^- Fortunately, the individual expressions 
for aE and gb are not needed to find a simple expression for the bulk viscosity. Returning 
to Eq. fl69ll we have 


c = 



dr: 


Pal feq/ 4par 
Ja \^a 




(93) 


Now the trick is to take a judicious linear combination of the conditions of £t. Add 
T — l^B times flS^ to This gives 


as 



+ dE 




dT 



The coefficient of qb is just {Oub/ dT)^, and from Eq. (I50ll the coefficient of qe is just 
{de/dT)cr- Therefore we have 



Since enters into Eq. (j9^ it is useful to derive the thermodynamic relations 


Tv 


2 

n 


W 

{de/dT)^ 


Ub 

{duB/dT)^ ■ 


(96) 


First, we derive the relation between the derivatives appearing in the above equations. Using 
de = Tds + jjiBduB and ds = UBda + aduB, we obtain 


de\ w 
Ub 


dT / 

/ o 


On 


B 


dT 


(97) 


Now for (duB/dT)^ we use Eq. (1^ . the third equality of Eq. (154)) . and the first equality 
of Eq. (15^ to obtain 




(98) 


Together with the previous equation we obtain the desired result fl96|) . Using these results 
in Eq. fl95l) we have 


Tn^as + TwaE = -vlV ^ J dr-J^Al‘X . (99) 

Making this substitution in Eq. flHHj) we obtain the expression ([92]). 

A similar arbitrariness arises in Eq. flH^ . Due to energy-momentum conservation, if we 
have a particular solution we can generate another solution as Ba = 7?^“ — b, where 
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6 is a constant independent of particle species a. This freedom is resolved by the Landau- 
Lifshitz condition of ht which requires that = 0 in the local rest frame. Starting with 
expression (IMj) we have 


= 




dTl^E, 


-(i?r-&)prA 


(f) 


re 
J a 


eq 


( 100 ) 


Factoring out the spatial derivative, and making use of the momentum space isotropy, we 
require that 




dr; 


' E* 


-E = 


dv: 


\v^ 

E* 


-EaBrr. 


eq 


( 101 ) 


The integral multiplying b is just 3Tw so that 


3Tw 

Substitution into expression fTfll) gives 



dr; 


b: 


,* 12 


e: 


-EaBrfa^- 


dr: 


I Pa 


E* 



ns Eg 

w 


sr/r ■ 


( 102 ) 


(103) 


There is no ambiguity in the solution to Eq. flH^ for Cg, so the expression for the shear 
viscosity fl^ is unchanged. 


VIII. RELAXATION-TIME APPROXIMATION 


At this point, it is convenient to derive the relaxation time approximation formulas 
for the shear and bulk viscosities and thermal conductivity. We start with the Boltzmann 
equation with the Chapman-Enskog expansion: 


dt 


= Cg 


(104) 


The left-hand side of Eq. 11041 is given by Eq. [78] while Cg can be found in Eq. [TS] In the 
energy-dependent relaxation time approximation 12|], we assume particle species a is out 
of equilibrium {(pg ^ 0) while all other particle species are in equilibrium (06 = 4>c = (j>d = 
0). Using Eq. [TS] the collision integral Cg greatly simplihes, and the Boltzmann equation 
becomes 


dt 


= Cg = 




Tg 


(105) 


where the relaxation time Ta{Ep) for species a is given by 


^a{El 


E 

bed 

E 

cd 

E 

be 


1 + 5, 


ab 


dTldT:dT:r,^W{a,b\c,d) 


dr* dr;^fr(a|c, d) 
dr* dTlf^^W{c\a,b). 


(106) 
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Next, we replace the left-hand side of Eq. 11051 using Eq. Into the right-hand side, 
we substitute (pa using Eq. |39l Then we equate terms on the left- and right-hand sides by 
matching tensor structures, and we obtain particular solutions for the functions Ba, and 


Ca from (pa: 


4par _ 
a 3^ 


Ip:p 


E* dT 


d f Ea Pc 


^par _ ' ^ 

a 


ha 


nsEa 


w 


^par _ a 

“ ~ 2TE* 


Finally, we substitute Eqs. I1071I109I into Eqs. 
relaxation time formulas: 


(107) 

, (108) 

(109) 

and 11031 and obtain the desired 


V = 


1 




ip:t 

Ef 


r.(E:)f’\ 


( 110 ) 




3 KubT 


MED 

' Ef 


ip:i"+ mxe:^ 


w 


d fE. 


l-^a 


f, 


eq 

a 5 


E / M 


IP, 


* 12 


E*2 


■TaiK 


nsEa 


w 


/, 


eq 

a ' 


(111) 


( 112 ) 


A few observations are in order. First, the transport coefficients computed with Eqs. 
1110111121 are strictly non-negative, as they must be. Second, this non-negativity is ensured 
by the squares in the integrands which came from enforcing the Landau-Lifshitz conditions 
of £t. (Recall the derivation of Eqs. I9^and fl03l i This shows that it is absolutely vital that 
the Landau-Lifshitz conditions are carefully enforced in order to obtain the correct results. 
A third point is that Eqs. IllOl and 11111 are obvious generalizations of the formulas obtained 
in previous works fl^. 271 to hnite baryon chemical potential. The crucial insight is that 
entropy per baryon [a = s/ub) is conserved in zeroth-order (ideal) hydrodynamics, so that 
variable that must be held hxed when deriving the variations from equilibrium. 


IX. CONCLUSION 

In this paper, we developed a flexible relativistic quasiparticle theory of transport co¬ 
efficients in hot and dense hadronic matter. A major goal was the simultaneous inclusion 
of temperature- and baryon chemical potential-dependent quasiparticle masses with scalar 
and vector mean helds, all in a thermodynamically self-consistent way. Classical statistics 
were used throughout to simplify the presentation, although complete results with quantum 
statistics are given in the appendix. From the dispersion relations for the quasiparticles, 
we derived the Boltzmann equation and then the transport coefficients using the Chapman- 
Enskog expansion. Next, we derived compact analytic expressions for the shear and bulk 
viscosities and thermal conductivity. These formulas can be used with the relaxation time 
approximation; alternatively, we have provided integral equations which may be solved for 
greater accuracy. We have shown that the transport coefficients are non-negative in the 

































20 


relaxation time approximation (as they must be) which is a direct consequence of carefully 
enforcing the Landau-Lifshitz conditions of ht. 

We also showed that previous bulk viscosity formulas (derived assuming zero baryon 
chemical potential) generalize straightforwardly to hnite baryon chemical potential if one 
recalls that entropy per baryon is conserved in ideal hydrodynamics. This was the crucial 
detail that allowed us to compute the variations from equilibrium and use them to derive 
the bulk viscosity and thermal conductivity formulas. 

It is a trivial matter to include a variety of scalar and vector helds, that is simply a matter 
of book-keeping. The same is true of additional conserved charges beyond baryon number. 
In future work we will study specihc hadronic models, including numerical solutions, along 
the lines of Ref. [l^ . 


Acknowledgment s 

This work was supported by the US Department of Energy (DOE) under Grant No. 
DE-FG02-87ER40328. 


Appendix 


This appendix has two goals. The hrst is to summarize the important results derived in 
the main body of the manuscript. The second is to include the effects of quantum statistics. 
All results presented here include quantum statistics. The limit of classical statistics is 
attained when |/a| <C 1. 

Departures from local kinetic and chemical equilibrium for particle species a are ex¬ 
pressed in terms of the function (pa as 




(113) 


We let 6 fa represent the deviation expressed in terms of the equilibrium energy while 6fa 
represents the deviation expressed in terms of the total nonequilibrium energy Ea] it is the 
latter which is conserved in local collisions and the one relevant for transport coefficients. 
The deviations to each other by 

5/. = 5/. + (1^) SE, = iA - (1 + dJD . (114) 

\(y^a J T0,IJ,% ^ 

Here the notation is da = (—1)^^“. We need to relate the variations in T and hb to the 
variation 6fa- The latter variation is done at hxed Ea and is 


Sfa 


f 

J ( 


eq 


77a ~ /J'a E T 



(1 + daD 


6T_ 

rp2 


(115) 


Here in what follows, the derivative is carried out at hxed entropy per baryon a. The factor 
from Eq. flll4p which needs to be rewritten is 


dEa 


fl^ (1 + dafl^) 



T{dEaldT)a 

^^a + T{^^^a/^T)^_ 


T 


(116) 
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In terms of 6fa the deviations in the energy-momentum tensor and baryon current are as 
follows. 


= 




5r« = /irf-’P-hA 

ad a 

(117) 

ad a 

(118) 

f, T(dEJdT), ) 

“ “ 1 Ea-fia + T (dfia/dT)^ / ' 

(119) 

sj‘B = Y^d 

(120) 

fdrh toe jar), i - 

J “1 Ea-fIa + T{dfia/dT)J 

(121) 


The collision term on the right side of the Boltzmann equation reads 

1 


"““Si s.. / 


dTldTldT*aW{a,b\c, d) 


X ^ fcfd (1 + dafa) (1 + — fafb (1 + dcfc) (1 + ddfd) 

+ J2f dTldTlW(a\cJ) 

cd 

x|/c/d (1 + dafa) — /a (1 + dcfc) (1 + ddfd) | 

+ ^ [ dTldTlW{c\a,b) 

be '' 

X |/c (1 + dafa) (1 + dbfb) — fafb (1 + dcfc) | • 


( 122 ) 


This expression explicitly includes 2 -H- 2 and 2 -H- 1 reactions. Higher order reactions are 
included in an obvious way. 

We now consider small departures from equilibrium, meaning that we keep terms only 
linear in the fa- We use chemical equilibrium; for example, a + b ■H- c + d gives 




/:vr ( 1 + daD ( 1 + dbO = /:vr (i+(i+ d^f, 

Then the collision integral becomes 

C„ = f dTldridr-,W(a,b\c,d) 

tt 1 + -I 

X {/r/r [(1 + d^fP) ^0 + (1 + rfo/;-) 4>4 - /r'/f [(i + d,fp) 4>. + (i + dji^) 4,,] 

cd d 

+ Y, j d:rldFlW(c\a,h)[- !^[(l + diC)<l>a + (l + dJT)4t] + (123) 
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The left-hand side of the Boltzmann equation is computed using the local equilibrium 
form of the distribution function. One form of the left-hand side (in the local rest frame) is 


dt 


= (1 + dafD 


+ (1 + daD 

- f:^ (1+ dj!^) 


IP, 


* 12 


3TE* 




dT 


dpUp 


“ w j E* ^ \T ) 


P^Pa 

2TE* 


DpUi/ -\- Ei/Up -\- ^^pi/dpuF 


(124) 


Now we subtract the right-hand side from the left-hand side and set the resulting ex¬ 
pression to zero. This leads to 

Aa (dpU^) + - cr (EpU, + E,Up + lAp^dpUf^) = 0 (125) 


where 
A = 


■\* 12 


3TE;! 




d f Ea fia 


T 


fr (1 + daf:^) 


E 

bed 


1 


1 + dab 


dTl dTl dT*aW{a,b\c, d) 


X {/r’/r [( 1 + difT) A.+( 1 +4/;’) A 4 

- f^fT 1(1 + dtf,") A. + (1 + i/D At] } 

+ dr:dniF(a|c,d){n(i+4/r)xio+(i+4r)4j-r/rA 

cd d 

+ Y^f dr:dr:iF(c|o,t){-4'’i(i+4/r')xi.+(i+4/r')xid + /:vr'xio}, 


(126) 


- 


nBEa\ 


w 


(1 + daf^) 


, . 1 + 
bed 


E 


dr*dr*dr*fT(a,6|c, d) 


X {/.“>/r 1(1+ djT) Btpi +(1+4r) bjpsi 
- /r/r 1(1 +dtD BtK +(1+40 firfi} 

+ Y. f dr;dnioic,4{/ri(i+4/7)B.p;i+(i+40B,ip'ji-/c'vrB.p^} 

cd d 

+ Y,f <ir;dr:i¥(c|o.6){-/f< 1(1 + 40 BaP^ + li + ^OBrfl + orBoPj}. 

be d 

( 127 ) 
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cr 


^flHl + dafD 

f drtdrtdr,w(,a,b\c,d) 


X {/:vr 1(1 + dJD CoP?p" + (1 + djp) 

- nvT [(1 + A/D CM + (1 + dM MM ] 

+ Y.f dTidr'M(Md){fl''[(i + dM)CM + (^ + dM)CM^-fT'fTCMM?j 

cd d 

+ Y.fdndr: W(c\a, 6){ - /;- [(1 + dM)CaPX + (1 + dj^)CM] + frJTCM} ■ 

be d 

( 128 ) 


Due to the tensorial structure of these equations the solution requires that = 0, = 0, 

and = 0. These are integral equations for the functions Aa, Ba, and Ca which depend 
on the magnitude of the momentum p*. 

The solutions for Aa and Ba are not unique. It is necessary to specify whether 
represents the flow of energy (Landau-Lifshitz) or baryon number (Eckart). We enforce the 
Landau-Lifshitz condition, sometimes known as the condition of £t, using any particular 
solutions. The results for the transport coefficients are 


C = 5E/< 




I Pa 


A = i 


W 


3 ynsT 

V = 


a 

2 

Is 


E / 

^ ^ a 

E / dr-j^cni 



4 par feq 

J a 5 

(129) 

a - —] Brfr‘ . 

(130) 

w / 



/Opar .req 

J a • 


(131) 


The particular solutions need not even satisfy the Boltzmann equation to satisfy the condi¬ 
tion of £t. 

A common approximation is the energy-dependent relaxation time approximation. It 
assumes that only one 0a is nonzero and the others vanish. Then the Boltzmann equation 
is approximated by 


dt 



(132) 


where the relaxation time Ta{E*) for species a is given by 


1 + dgf!^ 

rJE:) 
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bed 


1 + 


ab 


dr; dr: dr: w{a, b\c, d)fp (i + dj:^) (i + d.f. 


eq\ 


/ dr: dr, w (a|c, d) (1 + dj:^) (1 + d,f^ 

cd d 

E [ drldT:W{c\a,b)]f'{l + dM) . 




be 


(133) 
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The particular solutions are 
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* 12 
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E* 


nsEa 
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(1 + dafl^) , 
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“ “ 2TE: 


(1 + dafl^) . 


Substitution gives the transport coefficients 
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15T 
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Ef 


I *|2 I o 2^712 jTi* ^ (Ea Ho 

|Pal +3u„T E^ — 


1 2 
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3 XubT 

These are clearly positive-dehnite. 


E I‘‘E'J-^r^E:) (b. 


nsEa 
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C (1 + daC) ■ 
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